Draft version January 12, 2010 

Preprint typeset using I^T^^X style emulatcapj v. 11/10/09 



CONTRIBUTION TO THE EXTRAGALACTIC GAMMA-RAY BACKGROUND FROM THE CASCADES OF 

VERY-HIGH ENERGY GAMMA RAYS FROM BLAZARS 

TONIA M. Venters^'^ 
Department of Astronomy and Astrophysics, The University of Chicago, Chicago, IL 60637, USA 

Draft version January 12, 2010 



o 

(N 
oo 

o 



> 

m 
cn 



o 
o 



X 



ABSTRACT 

As very high energy (VHE) photons propagate through the extragalactic background light (EBL), 
they interact with the soft photons and initiate electromagnetic cascades of lower energy photons and 
electrons. The collective intensity of a cosmological population emitting at VHEs will be attenuated 
at the highest energies through interactions with the EBL and enhanced at lower energies by the re- 
sulting cascade. We calculate the cascade radiation created by VHE photons produced by blazars and 
investigate the effects of cascades on the collective intensity of blazars and the resulting effects on the 
extragalactic gamma-ray background. We find that cascade radiation greatly enhances the collective 
intensity from blazars at high energies before turning over due to attenuation. The prominence of the 
resulting features depends on the blazar gamma-ray luminosity function, spectral index distribution, 
and the model of the EBL. We additionally calculate the cascade radiation from the distinct spectral 
sub-populations of blazars, BL Lacertae (BL Lacs) objects and flat-spectrum radio quasars (FSRQs), 
finding that the collective intensity of BL Lacs is considerably more enhanced by cascade radiation 
than that of the FSRQs. Finally, we discuss the implications that this analysis and upcoming Fermi 
observations could have for the nature of the EBL, the evolution of blazars, blazar spectra, and other 
sources of gamma-ray radiation. 

Subject headings: diffuse radiation — galaxies: active — gamma rays: general 



1. INTRODUCTION 

The gamma-ray sky consists of resolved point sources 
(such as normal galaxies and active galaxies, pulsars, 
etc.), transient gamma-ray sources (e.g., gamma-ray 
bursts), and the diffuse gamma-ray radiation comprised 
of emission from the Galaxy and the extragalactic 
gamma-ray background (EGRB). The origins of the 
EGRB are, as yet, unknown; however, it is expected that 
the intrinsic emission of unresolved, extragalactic point 
sources comprises a sizable contribution to the EGRB. 
Additionally, many of these extragalactic point sources 
are also sources of very high energy (VHE) gamma ray^, 
which interact with the soft photons of the extragalac- 
tic background light (EBL) to create electromagnetic 
(EM) cascades, giving rise to another contribution to the 
EGRB. 

The gamma-ray sky, in general, was first observed by 
the SAS-2 and COS-B satellites in the 1970s. During 
the 1990s, the energetic gamma-ray experiment telescope 
(EGRET) aboard the Compton Gamma-ray Observatory 
observed the gamma-ray sky and, in particular, resolved 
271 gamma-ra y sources (the third c atalog of EGRET 
point sources; iHartmann etall[l999h of which 93 were 
identified, either confidently or potentially, as blazars 
(gamma-ray loud active galactic nuclei (AGNs)). Cur- 
rently, the Fermi Gamma-ray Space Telescope is taking 
data and in its first few months of observations has al- 
ready identified 109 blazars (LAT Bright AGN Sample 
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fLBAS): lAbdo et all [2009a[ ) . Thus, blazars constitute 
the class of gamma-ray emitters with the largest number 
of identified members. 

Blazars, specifically, are those AGNs for which the 
j et is closely a ligned with the observer's line of sight 
(jBlandford fc K onigl 1979), and include both BL Lac 
objects and gamma-ray-loud fiat spectrum radio quasars 
(FSRQs). Since resolved blazars comprise the class of 
gamma-ray emitters with the largest number of identi- 
fied members, unresolved blazars are expected to have 
a siz able contribution to the EGRB (i Sreekumar et al.1 
|1998( ). However, the exact amount of this contribu- 
tion is, as yet, unclear due to the uncertainty sur- 
rounding the distribution of blazars with respect to lu- 
minosity and redsh ift, the blazar gamm a -ray luminos- 
ity fu n ction (CLE: iPadovani et ah Il993t IStecker et al.l 



niga'; 'Salamon fc Stecker' ' 19941: Ich iang et al.' '19951: 



.Stecker fc Salamon 1996; Kazanas &: Pe rlman 1997 
Chiang fc M ukhcriec 1998; Mukheri ee fc Chiand 11999 : 



Mucke&Pohl 2000: .Dermer. .2007: Giommi et al.l 



2006 



Narumoto fc Totan^ l2006l hereafter NT06). 



In addition to the dependence on the blazar GLF, the 
blazar contribution to the EGRB also depends on the 
nature of blazar spectra at GeV energies. The energy 
spectrum of the gamma-ray emission of blazars at these 
energies has traditionally been approximated by a power 
law, defined by a spectral index, which, in most cases, has 
a value between 1.5 and 3 for individual blazars. The dis- 
tribution of blazar gamma-ray spectral indices, the spec- 
tral index distribution (SID), is a critical input in the de- 
termination of the shape of the spectrum of the collective 
unresolved blazar i ntensity t hat is the blazar contribu- 
tion t o the EGRB (IStecker fc Salam ol]ll996l: iPohl et al.l 
[19971: iPavhdou fc Venters! 12008c i Venters et al.l I2009D . 
The spread in the blazar SID determines the fraction 
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of blazars with hard spectra, which wiU contribute most 
significantly at high energies and hence will introduce 
curv ature i n the shape of the unresolved blazar inten- 
sitv dStecker & Salam onHTOOl FVe ntcrs & Pavhdou 2007; 
iPavlidoufcVenterd 120081 : 1 Venters et al.l i2009). Since it 
is the highest energy photons that interact with EBL 
photons, the determination of the blazar SID is crucial 
to our understanding of the contribution of EM cascade 
radiation to the EGRB. 

The EBL is composed of emission from starlight (at 
optical, ultraviolet, and near-infrared wavelengths) and 
reradiated thermal dust emission (at far-infrared) in 
galaxies. At observed energies beyond the EGRET 
energy range (but well within the Fermi energy 
range), photons suffer significant attenuation due to 
pair produ ction interactions with t h e soft photons of 
the EB L (Salamon & Stcckcr' 1998; Chen ct al. '20p; 
Kneiskc ct al. 2002, 2004: Kneisk e 2008: Stecker et al.i 
2006', 2007; Prim ackeTall |200S [Gilmoreet all 120091: 



Venters et al.ii2009i ). In lVenters et all (|2009f ). wedemon- 



strated that blazars exhibit an absorption feature at the 
highest energies in their collective spectrum. We showed 
that the strength of such an absorption feature depends 
on the blazar GLF with the strongest feature arising for 
those models that situate more high-luminosity sources 
at high redshifts. We also showed that the shape of the 
absorption feature depends on the EBL model. Thus, 
we demonstrated that if blazars dominate the EGRB, 
the measurement of the absorption feature at the highest 
energies can place constraints on the blazar GLF and the 
EBL. We did not, however, propagate high energy pho- 
tons to determine the resulting cascade radiation, which 
will further alter the shape of the collective unresolved 
blazar spectrum. In this paper, we study the effect of 
the cascades resulting from propagation on the collective 
blazar intensity. 

In interacting with EBL photons, VHE gamma rays 
will produce pairs of electrons and positrons, which will 
inverse Compton scatter EBL photons to high energies. 
These upscattered photons will, in turn, interact with 
soft EBL photons through pair production, and this "EM 
cascade" process continues until the energies of the re- 
sulting photons are low enough that pair production is no 
longer efficient. For any cosmological population emit- 
ting gamma rays at VHEs, the effect of EM cascading 
results in a flux suppression at the highest energies and 
enhancement at lower energies. In the case of blazars, 
which can comprise a sizable contribution to the EGRB, 
predictions of the resulting enhancement at lower ener- 
gies can lead to the overproduction of the EGRB if the 
collective high-energy intensity of blazars is high and/or 
the EBL is high (Coopi & Aharonian 1997). M ore re- 
cently , iKneiske fc MannheimI ([20081) and llnoue fc Totanil 
()2009f ) have estimated the contribution of cascade radi- 
ation to the EGRB. They insightfully demonstrate that 
in including radiation from EM cascades, blazars can ac- 
count for nearly all of the EGRB (- 80%-90%). However, 
in neither analysis was the blazar SID as determined from 
either EGRET or Fermi data included, and thus, the im- 
pact of the population of blazar GeV spectral indices was 



not fully examined. Furthermore, neither analysis makes 
use of a Monte Carlo propagation scheme, and hence do 
not fairly sample secondary electron energies. 

In calculating the spectrum of cascade photons, we 
make use of a Monte Carlo propagation code called Cas- 
cata, which samples secondary electron energies from 
cross-section- weighted distributions (see Section 12.21 and 
the Appendix). Once sampled, the resulting spectrum of 
EM cascade radiation is then calculated. In this manner, 
Cascata calculates a more fully sampled cascade spec- 
trum resulting from a spectrum of primary photons. 

In this paper, we revisit the contribution of EM cascade 
radiation to the EGRB. Specifically, we study the impact 
of the blazar GLF, the blazar SID, and the EBL model 
on the spectrum of cascade radiation, which, in turn, af- 
fects the blazar contribution to the EGRB. All of the 
aforementioned inputs remain quite uncertain. However, 
if blazars do, in fact, comprise the bulk of the EGRB in- 
tensity, then the careful study of EM cascades along with 
Fermi observations can be used to constrain the inputs 
of the collective intensity and the nature of the EBL. 
Furthermore, as the EM cascade radiation is sensitive to 
the blazar SID, it is also sensitive to possible breaks in 
blazar spectra beyond the EGRET energy range (e.g., 
spectral breaks or cutoffs). Thus, the study of EM cas- 
cades together with observations of the EGRB can also 
provide information about blazar spectra beyond tens of 
GeV, which can, in turn, provide insight into the nature 
of blazar emission. 

In this paper, we demonstrate that with a Monte Carlo 
propagation code such as Cascata and independent con- 
straints on the inputs from Fermi observations, as well 
as those from imaging atmospheric Cherenkov telescopes 
(lACTs), the blazar contribution to the EGRB wih ul- 
timately be determined. In Section 2, we present the 
formalism of the calculation of the collective unresolved 
blazar intensity and a short discussion of the aspects of 
the code relevant to the propagation of VHE gamma 
rays. In Section 3, we discuss the inputs of the calcu- 
lation and their uncertainties. In Section 4, we present 
the results of the calculation, and we discuss these results 
in Section 5. 

2. FORMALISM 

2.1. Contribution to the ECRB from Cascades of VHE 
Photons from Blazars 

The contribution to the EGRB due to blazars can be 
viewed as the superposition of the collective intensity of 
intrinsic blazar spectra and the contribution from the 
cascade radiation from the interactions of VHE photons 
from blazars with the EBL: 



Ie{Eo) = pEiEo) + I%{Eo 



(1) 



where the intensity, Ie{Eq)^ is given in units of pho- 
tons per unit energy per unit time per unit area per unit 
solid angle emitted at observer frame energy, Eq. The 
collective intensity of intrinsic blazar spectra including 
attenuation by the EBL is given by (for derivation, see 
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where Ef is some fiducial observer 
frame energy (taken to be 100 MeV), 
= EfdN{E'j)/dtdE is the isotropic gamma-ray 
luminosity of a blazar at the fiducial rest frame energy, 
E'f = Ef{l + z'), = d'^N/dL^dVcom is the blazar GLF, 
t(Eq,z) is the optical depth as a function of observed 
photon energy and source redshift (see Section 13. 3p . 
Pi (a) = dN/da is the SID of gamma-ray blazars cor- 
rected for measurement error and sample bias (see Sec- 
tionESD, and L^,„,ax = 47rd|(a - 1)(1 + z)°'-^EfF^^^in, 
where i^-y,min is the sensitivity of the survey above the 
fiducial observer frame energjQ. In deriving the above 



equation, we have assumed the standard ACDM cos- 
mology and that blazar energy spectra can be described 
as single power laws defined by the spectral index, a 
{dN^/dE^ oc £;-"). 

Making the dependencies explicit, the cascade intensity 
is given by 



d'^N'= 



dtdAdndE ' 



(3) 



where dN'^/dE is the spectrum of cascade photons due 
to pair production and Inverse Compton scattering (for 
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where dNr/dE is the spectrum of inverse Compton 
scattered radiation per electron of Lorentz factor, F, 
P(/; Ep, z) is the probability that the pair production 
interaction at a given redshift, 2;, of a primary photon of 
energy Ep will produce electron-type particles of energies 
Eel = fx Ep and Ee2 = {1 - f) x Ep, Ti = f x Ep/mc^, 



where A'zbins is the number of redshift bins, zi is the red- 
shift of the l^^ bin, Azi is the size of the l^^ redshift bin, 
N^{zi) is the number of photons at zi undergoing pair 
production, Nj^ = 2 is the number or electrons produced 

by fc"^ photon undergoing pair production, and dA^p*; /dE 

j 

is the spectrum of inverse Compton photons produced by 



r2 = (1 - /) X Ep/mc , and d N^/dzdEp is the con- ^j^e j"^ electron with Lorentz factor T*^. We appropriately 
tinuous spectrum of photons undergoing pair production 
interactions as a function of redshift and primary energy. 

While numerical integration of Equation ^ is, in prin- 
ciple, possible, doing so would be both complicated and 
computationally expensiv^^. Hence, it is desirable to 
make use of a Monte Carlo method, such as Cascata, 
that propagates photons from known spectra and red- 
shift distributions and calculates the cascade spectrum 
of each photon, resulting in an observed spectrum char- 
acterized by a discretized version of Equation (j4]) (for 
derivation, see the Appendix), 



normalize the resulting spectrum in the final step. 

In determining the cascade contribution to the EGRB 
from VHE photons from blazars, primary photon ener- 
gies are randomly generated according to the absorbed 
EGRB intensity (in observer's quantities), which is the 
differential amount of radiation that has been "absorbed" 
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This intensity is interpreted as the radiation that creates 
EM cascades. As such, it will also be used to normalize 
the cascade spectra in the final step. 

For EGRET, the sensitivity above 100 MeV is lO"'^ 
pliptons cm~^s~^. For Fermi LAT, the one-year sensitivity above 
Furthermore, Equation l(4]| is only valid for the case of EM 
cascades developed through pair production and inverse Compton 
interactions. While similar equations could possibly be derived for 



The absorbed EGRB intensity is also used to gener- 
ate the redshifts of the sources of the VHE photons by 



100 MeV is 2 X lO"'' photons cm- 
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The created photons are then propagated until the red- 
shift of interaction, randomly generated from the proba- 
bility of interaction, p = 1 — exp[— r(iJo, 2^)]- 

2.2. The Cascade Formalism of Cascata 

In addition to generating primary photons, background 
photons with which the primary photons interact must 
also be gene rated according to the appropriate distri- 
butions (see iProthero"^ 119861 for details) . For a given 
redshift, z, the background photon energies, e, may be 
sampled according to 



p(e) = Xj{Ep, z)n{e, z) 



(1 — coa9)d{cos9) , 



(8) 

where X-y{Ep^z) is the mean free path of a gamma- 
ray photon of energy, Ep, n(e, z)de is the proper pho- 
ton number density (EBL plus cosmic microwave back- 
ground (CMB)), 9 is the angle between the paths of 
the gamma ray and the background photon, and c^-yis) 
is the pair production cross section as a function of 
the square of the center of momentum (CM) energy, 
s = E'^j^ ~ 2eEp{l — cos 6). The mean free path (see 
Figure [T]0, \~f{Ep, z), of the gamma-ray photon is given 
by 



\-z\Ep,z) = 



^ [l-cosl 



n{e, z)aj^{s)ded{cos 9) . 

(9) 

where eth = 2m1c^ / Ep{l — cos9). The pair production 
cross-section is given by 
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(10) 

where (3c = (1 ~ ^m'lc'^ / sY^'^c is the magnitude of the 
velocity of the electrons in the CM frame. The angle be- 
tween the paths of the two photons can be determined by 
sampling the square of the CM ener gy, s, from Equation 
([U differentiated with respect to s (|Protheroe.l98 6f1 



(11) 



The scattering angle, a, of the secondary electrons in 
the CM frame is determined from the differential pair 
production cross section: 
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(12) 



® The attenuation for photons of energies less than 10^^ eV is 
dominated by the EBL. The various EBL models predict different 
amounts of UV radiation resulting in the differences seen in the 
figure in the mean free paths at energies of interest for lACTs and 
Fermi. 

In differentiating Equation Jsjl with respect to s, we are, in 
fact, differentiating only with respect to cos 8, where again s = 
2eEp{l - cos 9). 
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Fig. 1. — Pair production mean free path for VHE photons 
propagating through the present-day (z = 0) soft photon back- 
ground as a function of energy for four different models of the 
background. Solid: Soft photon background consisting of only the 
CMB. Dotted: Soft photon background consisting of the CMB and 
the Gilmorc ct al. 2009 model of the EBL. Dot-dashed: Soft pho- 
ton background consisting of the CMB and the Stcckcr ct al. 2003 
model of the EBL. Dashed: For comparison, the soft photon back- 
ground consisting of the CMB and the iKneiskai2008i model of the 
EBL. 



where C = cos a. The energies of the electrons in the 
lab frame are then completely specified by the energy 
of the propagating gamma ray, the randomly generated 
parameters of the pair production interaction (e, 9, and 
a), and relativistic kinematic^ 



and 



Eel = 7C 



Ee2 = 7C 
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+ cos(a -I- ip) 



■/3|Pe2 1 cos(a + ■(/') 



(13) 



(14) 



where 7 and (3 determine the boost from the lab frame 
to the CM frame, ip is the angle of the photons with re- 
spect to the X-axis in the CM frame, = [E^^/Ac"^ ~ 



mgC^)^/^ and jpgal = {E^yil^c? -mlc^)^/'^ are the mag- 
nitudes of the CM momenta of the secondary particles. 

Finally, the electrons scatter background photons 
through the inverse Compton process resulting in a scat- 
tered photon spectrum given by (see the Appendix) 



dNr 

'dW 



d^Nr' 
dtdEde 



dt 



dE, 



dEf.de , 



(15) 



where d'^N^' /dtdEde is dif feren t ial scattered photon 
spectrum given by Equation (jA2|) ([Blumenthal fc Gouldl 



° The procedure for determining the energies of the two sec- 
ondary electrons created through the interaction between the two 
primary photons involves boosting from the lab frame to the CM 
frame, determining the energies and momenta of the two secondary 
electrons in the CM frame, and then boosting back to the lab frame. 
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119701) and dEe/dt is the electron loss rat43- 
3. INPUTS 

3.1. Blazar Gamma-ray Luminosity Functions 

Determining the blazar GLF from observations re- 
lies on the ability to associate gamma-ray blazars with 
lower energy counterparts for which redshifts can be mea- 
sured. However, identifying gamma-ray-loud blazars can 
be complicated due to the large positional error circles 
of GeV gamma-ray experiments. Blazar variability and 
low number statistics further complicate the determina- 
tion of the GLF. Thus, while the blazar GLF is probably 
one of the most studied and debated properties of the 
gamma-ray-loud blazar population, it remains, to this 
day, uncertain. 

In gamma-ray observations, blazar GLFs are con- 
structed from luminosity functions based on survey data 
taken at wavelengths for which at least some of the obser- 
vational uncertainties are not as challenging. The most 
uncertain parameters (e.g., normalization due to rela- 
tivistic beaming and the faint- end slope) are then fitted 
to gamma-ray data (see e.g., IStecker fc SalamonI 119961 : 
IGiommi et al.l [20061 NT06). Such an approach raises a 
number of concerns: 

(1) Not all blazars are necessarily explained by the 
same emission process, and different types of 
blazars (i.e., BL Lac- like objects and FSRQs) could 
form separate populations with res pect to emission 
(jSikora et al.l 120021: iBottcheii [20071) and, hen ce, re- 
quire separate lumi nosity functions ((Muckc fc Pohll 
[2000l:rDermed[2Q07l) . 

(2) Flaring blazars (and different types of flaring 
blazars) could also form separate populations with 
respect to emission and require sepa rate luminosity 
functions (jStecker fc Salamonlll996D . 

(3) Many unidentified gamma-ray sources could be 
unidentified blazars, which, without identification, 
would result in underestimating the normalization 
of gamma-ray blazars with respect to low-energy 
blazars and the introduction of uncertainties in the 
redshift distribution of resolved blazars. 

Alternative approaches have been proposed (see, e.g., 
[Permer 2007.) . but they are also subject to uncertainty. 

With the availability of Fermi data, many more blazars 
will be observed, and at least some of the aforementioned 
uncertainty will be alleviated. In fact, with the release 
of the first year data approaching, predictions for the 
Fermi source counts will be directly tested. Further- 
more, if the collective intensity of blazars does, in fact, 
constitute a substantial contribution to the EGRB, then 
Fermi will o bserve an absorptio n feature at the high- 
est energies (j Venters et al.l I2009D . The shape of such 
an absorption feature is highly sensitive to the blazar 
GLF and, hence, can constrain the evolution of blazarif^. 

^ Note that this process must be repeated iteratively if VHE 
electrons are produced. 

It should be noted that the shape of the absorption feature 
is also highly sensitive to the EBL (see Section |3.3|I through which 
VHE gamma rays propagate. As such, the absorption feature can 
also constrain EBL models. 



While Fermi data from the first fe w months of o bserva- 
tions of AGNs have been released ()Abdo et al.| [2b09a). 
the sources listed in the catalog represent high sensitiv- 
ity detections (~ lOcr) and do not constitute a complete 
sample. Thus, for the purposes of this paper, we make 
use of the best-fit pure luminosity evolution (PLE) and 
luminosity-dependent density evolution (LDDE) models 
of the likelihood analysis of NT06. Calculations based 
on luminosity functions, based on Fermi first year data, 
will be addressed in a future publication. 

3.2. Blazar Spectra 

As with the GLF, determining the contribution of 
the collective intensity of any population of gamma- 
ray emitters to the EGRB requires precise knowledge 
of their spectra. If all blazar spectra were described 
by the same power law, then the unabsorbed collec- 
tive intensity would also follow that same power law. 
If, on the other hand, the spectral indices of the pop- 
ulation of blazars form a distribution with a given fi- 
nite spread (the SID), the spectrum of the unabsorbed 
collective intensity will exhibit curvature due to the in- 
creasing relative importance at high energies of blazars 
with harder spectral ind ices (jStecker fc SalamonI 119961 : 
IPavlidou fc Venter^ '2(Wl. Since the radiation at the en- 
ergies which are likely to initiate EM cascades is most 
sensitive to the effect of the spread in the SID, it is es- 
sential that the blazar SID be determined. 

In the analysis of blazar spectral indices, we make 
the critical assumption that the spectral index of a 
blazar is independent of its luminosity (as indicated in 
I Venters fc Pavlidou[ 120071 ). Under this assumption, the 
incompleteness of the Fermi LBAS with respect to lu- 
minosity does not hinder the determination of the SID 
even though it complicates the extraction of blazar GLFs. 
Thus, for the purposes of this paper, we determine the 
blazar SIDs from the Fermi LBAS blazars foll owing 
the lik elihood procedure outlined in Venters fc Pavlidod 
(|2007l ). For the Fermi LBAS blazars, we find the pa- 
rameters of the maximum likelihood (Gaussian) SID to 
be (ao — 2.23, co = 0.26). For the BL Lacs, we find 
(ao = 1.97, do = 0.175); and for the FSRQs, we find 
(ao — 2.37, (Jo — 0.14j3- We also correct for the sam- 
ple bias introd uced in using a flux -limited catalo^3 (for 
discussion, see lVenters et al.|[2009() . 

Finally, while it may be true that blazar spectra ap- 
pear as single power laws in the GeV band, it may not be 
so beyond the GeV band (or, at least, it may not be the 
same power law) . Observations by lACTs demonstrate 
that blazar spectra are steeper at VHEs ([Albert et al.l 
[200l) . and even the Fermi LBAS provides a few examples 

Note th at the best-f it sprea ds are slightly small er than that 
indicated by lAbdo et al.l I l2009al ). This is because lAbdo et al.l 
( 2009a) calculate the rms of the observed blazar spectral indices, 
which does not account for measurement uncertainty in each indi- 
vidual spectral index. 

Blazars observed in a flux-limited survey will be biased to- 
ward high-luminosity, low-redshift blazars. Low-luminosity, high- 
redshift blazars that are most likely to appear in a gamma-ray 
survey are those with spectral indices that are hard with respect 
to the population as a whole. However, in this case, the correc- 
tions will only be approximate since we are not calculating them 
from luminosity functions based on LBAS blazars. Furthermore, 
the LBAS is not exactly a flux-limited survey due to the variatio n 
in the diffuse intensity with galactic latitude IIAbdo et al.ll2009bj ). 
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of blazar spectra that appear to break at energies greater 
than ~ 10 GeV. In the calculation of the blazar contribu- 
tion to the EGRB including cascades, one would ideally 
include the effects of possible breaks in blazar spectra. 
However, doing so is quite complicated as it is not clear 
whether the breaks are due to changes in the spectra 
endemic to the sources or, at least in part, to attenua- 
tion, particularly for high-redshift blazars. Even if the 
breaks do reflect intrinsic changes in the source spectra, 
the cause(s) of the breaks is not, as yet, well understood, 
and hence the intrinsic spectral breaks and the energies 
of these breaks are uncertain. Additionally, comparisons 
between blazar GeV spectral indices and TeV spectral 
indices require simultaneous observations of blazars by 
Fermi and lACTs, and as yet, the sample size for such 
observations is quite smalO- Finally, the observations 
provided by the lACTs and the spectral breaks demon- 
strated in the LBAS do not constitute a flux-limited cat- 
alog, and as such, the implications of these observations 
for the blazar population as a whole remain unclear. In 
considering spectral indices from a flux-limited catalog, 
correcting for the observational bias introduced by the 
overrepresentation of hard er spectral indi ces in the SID is 
relatively straightforward (jVenters et al.l 12009). but this 
is not the case with a catalog that is not flux-limited. 
Thus, in light of the uncertainty surrounding blazar spec- 
tral breaks, we consider the study of their effects to be 
out of the scope of this paper; hence, for the purposes 
of this analysis, we assume that blazar spectra remain 
unbroken at VHEs. As such, the calculated cascades will 
be larger than expected for blazar spectra which soften 
at VHEs. 

We will return to the question of blazar spectral breaks 
in an upcoming publication. In the meantime, we note 
that Fermi measurements of the EGRB and theoretical 
determinations of the blazar contribution to the EGRB 
that include cascade radiation can provide constraints on 
blazar spectra at VHEs within the context of the blazar 
population as a whole, which would be, in effect, compli- 
mentary to the constraints provided by lACTs on indi- 
vidual blazars. 

3.3. Extragalactic Background Light 

The spectral energy distribution of the cosmic back- 
ground of infrared, optical, and ultraviolet light (the 
EBL) encodes much information about the astrophys- 
ical and cosmological processes in the universe. The 
EBL intensity originates from starlight (at optical, ul- 
traviolet, and near-infrared wavelengths) and reradiated 
thermal dust emission (far-infrared) in galaxies. Hence, 
the EBL is directly connected to the evolution with cos- 
mic time of structure in the universe. However, the na- 
ture of key contributors to the EBL is largely uncertain 
(e.g., the evolution of the cosmic star formation rate; 



the amount of UV radiation from young, massive stars 
that escapes from the surrounding gas; the amount of 
re-emission by dust in galaxies). Furthermore, observa- 
tions are hindered b y emission from our ow n galaxy (for 
detailed review, see iHauser fc Dwg[200l . Neverthe- 
less, despite the degree of complication, there are several 
approaches to overcoming the observational and theoret- 
ical uncertainties and, ultima tely, modeling the E BL (for 
more detailed discussion, see iKneiske et all [2001 . 

(1) Forward evolution of structure for mation and the 
resul t ing luminosity funct ions (e.g.. iPrimack et all 
[20081 : TGilmore et al1[2009h : 

(2) Backward evolution of observed local luminos- 
ity functions to higher reds hifts, most notably 
at infrared w avelengths ( e.g., [S teckcr e t all 120061 : 
IKneiske et all 12002,., 200l [KnSskc .2008) ; 

(3) Evolution determined from direct observations 
over some range in redshift such as the global 
star formation rate inferred from deep galaxy 
surveys or tracers of cosmic chemical evolution 
convolved with galaxy spectral energy distribu- 
tions determi ned from population synthesis models 
fe^Stecker et al.ll2006l: IKneiske et alMO^ [200l 
IKneiskdl2008fr 

While progress is being made on all fronts, each 
method remains, as yet, subject to its own uncertain- 
ties, so it is not obvious which is the best one to use in 
propagation studies. Thus, we adopt two models that 
are representative of the cur rent state of the fiel d and 
bracket the possibilities : the [ Stecker et al.l (|2006l ) EBL 
model and the Gilmor e et al.l (I2009D E BL model (for 
more information, see Vent ers et al.ll2069l and references 
therein). The most notable distinction between these 
two models is the amount of UV radiation predicted. 
Since the Stecker model does not include the effects of 
dust extinction, it predicts a higher amount of UV radi- 
ation than the Gilmore model. Other pr evalent models 
in the liter ature, such the [Kneiskd ( !2008 [) high-UV EBL 
model, the IKneiske (I2008D best-fit EBL model, and the 
iFranceschini et al.l ( 2008| rEBL model, yield results sim- 
ilar to those of the Stecker and Gilmore models (at least 
in the analysis presented in this paper). 

In propagating VHE gamma rays, we make extensive 
use of the optical depth of VHE photons due to pair pro- 
duction interactions computed from the proper photon 
number density, n(e, z)de: 

t{Eo,z)= dz'— sm9d9 de n {e, z') a{s) , (16) 
Jo dz'Jo J,^^ 



(1 + z'Y (1 + VLmz') -z'{2 + z') 



-1/2 



(17) 



fdz'^^rdz' ; , 

h dz' 7o Ho{l + z') 

However, there is considerable interest within the gamma-ray 
community to perform such coordinated observations, and even 
though the different natures of the involved experiments ensure 



that they will be performed on a case-by-case basis, the sample 
size of multi-wavelength observations will grow with time. 
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is the look-back time. 

4. RESULTS 

In this section, we discuss the effects of EM cascade 
radiation from the propagation of VHE photons from 
blazars through the EBL on the overaU coUective blazar 
intensity and the imphcations for their contribution to 
the EGRB. In each plot, we show the overall collective 
blazar intensity, the contribution to the overall collective 
intensity from the collective intensity of intrinsic blazar 
spectra, and the contribution to the overall collective in- 
tensity from cascade radiation arising from propagation. 
We also discuss the dependence of the cascade radiation 
on the blazar GLF, the blazar SID, and the model of the 
EBL. 

The blazar contributions to the EGRB including cas- 
cade radiation for two models of the blazar GLF are plot- 
ted in Figures [5] and |31 For both GLFs, the calculation 
was performed for two separate models of the EBL (solid 
lines represent the total intensity from blazars, including 
cascade radiation, for each of the four possibilities). Fig- 
ure [2] presents the calculation assuming the NT06 PLE 
model of the GLF, and Figure |3] presents the calcula- 
tion assuming the NT06 LDDE model of the GLFQ. 
The top pan e ls pre sent the calculations assuming the 
iStecker et al.l ()2006[ ) EBL model, while th e bottom pan- 
els present the calculations assuming the iGilmore et al.l 
(p009) EBL model. All calculations are made assuming 
the parameters of the SID calculated based on the Fermi 
LBAS (see Section IX^. Fo r com parison, we have also 
plotted the lSreekumar et al.l ()1998l ) determination of the 
EGRB as measured by the EG RET (square d ata points 
with statistical error bars), the lStrong et al.l (2004) de- 
termination of the EGRET EGRB (circular data points 
with statistical error bars), a gray band representing 
the Fermi preliminary EGRB (including statistical and 
instrumental systematic uncertainties) from the first 9 
months of dateP^. the collective intensity of intrinsic 
blazar spectra assuming no absorption (dotted lines), the 
collective intensity of intrinsic blazar spectra including 
only absorption of the highest energy photons (dashed 
lines), and the cascade radiation (dot-dashed lines). In 
this analysis, we include intrinsic emission up to 100 TeV. 

As evident, the amount of cascade radiation created 
in the propagation of VHE photons can be substantial, 
depending on the shape of the collective intensity of in- 
trinsic blazar spectra (for examples of cascade radiation 
resulting from other shapes, see below). For the shape 
determined by the SID indicated by Fermi data, the cas- 
cade radiation significantly alters both the amount and 

^'^ Note that for the best-fit parameters of these luminosity func- 
tions the collective intensity of blazars (not including absorption 
and cascade radiation) comprises ~ 50% of the EGRB. However, in 
several of the cases presented in NT06, parameters for which unre- 
solved blazars can account for 100% of the background are within 
the Ic contours. It should further be noted that with the parame- 
ters provided by NT06, the LDDE model predicts an overall lower 
intensity than the PLE model. 

There is some uncertainty in the determination of the param- 
eters of the blazar spectral index, which will result in uncertainty 
in the overall shape of the spectrum of the unresolved blazar inten- 
sity. The already greatly improved blazar number statistics in the 
Fermi LBAS has reduced this uncertainty; hence, for the purposes 
of this paper, we simply analyze the best-guess spectra. 

As presented by M. Ackermann at TeV Particle Astrophysics 

2009 
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Fig. 2. — The blazar contribution to the EGRB including EM 
cascade radiation from the propagation of VHE photons through 
the EBL assuming the PLE model of the blazar GLF. Top: The 
calculation assuming the Stecker EBL model. Bottom: The cal- 
culation assuming the Gilmore EBL model. Solid lines: The to- 
tal blazar contribution, including EM cascade radiation. Dashed 
lines: The collective intensity of intrinsic blazar spectra without 
cascades but including absorption. Dot-dashed lines: The inten- 
sity from EM cascades, including absorption. Dotted lines: The 
collective intensity not including absorption and cascade radiation. 
Square data points: Srcekumar ct al. (1998) EGRET EGRB. Cir- 
cular data points: ,Strong et al.. (2004 ) EGRET EGRB. Grey band: 
Preliminary Fermi EGRB with statistical and instrumental sys- 
tematic uncertainties. 



the spectral shape of the overall intensity. In fact, for 
model GLFs that predict the contribution of intrinsic 
blazar spectra to comprise only ~ 50% of the EGRB, the 
addition of cascade radiation greatly increases the overall 
contribution from blazars, at least at the highest ener- 
gies. At lower energies, the cascade radiation drops off 
slowly over several decades in energy, and the collective 
intensity of intrinsic blazar spectra rises quickly; thus, 
the contribution of the cascade radiation to the overall 
collective blazar intensity is much less at lower energies. 

Also seen in the figure is the contrast in appear- 
ance between the PLE model of blazars and the LDDE 
m odel, most notably a t the highest energies. As noted 
in lVenters et al.l ()2009( ). the collective blazar intensity is 
much more distributed with respect to redshift in the 
PLE model than in the LDDE model, in which the ra- 
diation tends to be concentrated at lower redshifts (see 
their Figure 2). Hence, just as the absorption feature in 
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Fig. 3. — Same as Figure [2l except for the LDDE model of the 
blazar GLF. 



the collective intensity of intrinsic blazar spectra is much 
more prominent in the PLE model, so is it in the cascade 
radiation calculated from this model. Thus, not surpris- 
ingly, while there is more cascade radiation below 10 GeV 
in the PLE model than in the LDDE model, above 10 
GeV, the cascade radiation in the LDDE model overtakes 
that of the PLE model. Since there is more high-energy 
absorption of the collective intensity from intrinsic blazar 
spectra in the PLE model, the cascade radiation in the 
PLE model dominates the intensity at energies beyond 
the EGRET energy range, contributing nearly 4 times as 
much radiation as the contribution from intrinsic blazar 
spectra at 100 GeV. Not so for the LDDE model, in which 
the cascade radiation contributes at most 50% more than 
that of the intrinsic blazar spectra at 100 GeV. 

In addition to considering various blazar GLF mod- 
els, we also investigate the impact of the EBL model. 
In both models, the effect of the cascade radiation is to 
flatten the overall collective intensity from blazars. How- 
ever, at the highest energies, there is less absorption in 
the Gilmore model than in the Stecker model. Below 10 
GeV, the cascade radiation for the two models is roughly 
comparable, though there is more radiation at the lowest 
energies for the Stecker model. Above 10 GeV in the PLE 
model of blazars, there is more cascade radiation for the 
Gilmore model than for the Stecker model due to the en- 
hanced absorption in the Stecker model. As a result, the 
spectra at high energies are vastly different. As noted in 



iVenters et al.l ()2009D , such an absorption feature is much 
more prominent in the Stecker case than in the Gilmore 
case owing to the higher UV background in the Stecker 
EBL model. Since the pair production cross section as 
a function of the CM energy peaks at twice the electron 
mass, one would expect that gamma-ray photons of en- 
ergies ^ tens of GeV are most likely to interact with UV 
background photons. Thus, unsurprisingly, models with 
high UV backgrounds will result in more suppression at 
high energies. On the other hand, in the LDDE model, 
while there is still more cascade radiation at high ener- 
gies for the Gilmore model than for the Stecker model, 
the resulting distinction between the two models is less 
prominent than in the PLE model owing to the prefer- 
ential distribution of relatively high-flux blazars at low 
redshifts. 

We also investigate the possibility that BL Lacs and 
FSRQs form separate populations with respect to emis- 
sion and evolution by propagating photons from distinct 
GLFs and SIDs for the two populations (Figure |4]). As 
both the third EGRET catalog and the LBAS indicate 
that the BL Lacs tend to be more concentrated at low 
redshifts as compared to the FSRQs, we assume their 
evolution to follow the LDDE model. By contrast, we 
assume the FSRQs to follow the PLE modeO For the 
SIDs, we assume the parameters determined from the 
Fermi LBAS for the respective sub-populations (see Sec- 
tion [XI]). 

The hard spectra of BL Lacs cause the resulting cas- 
cade radiation to be considerable relative to that of the 
collective intensity from intrinsic spectra even though 
their collective radiation tends to be concentrated at low 
redshifts in this model. In fact, the total intensity from 
BL Lacs could dominate the EGRB at high energies and 
even overproduce it, though it does not explain the lower 
energy intensity. By contrast, the relatively steep spectra 
of FSRQs are not conducive to much cascade radiation 
even though they are preferentially situated at high red- 
shifts, and most of the intensity attributed to FSRQs oc- 
curs at lower energies. Hence, while the BL Lacs (includ- 
ing cascades) dominate the collective blazar intensity at 
the highest energies, the FSRQs dominate at the lowest 
energies. In Figure [5l we have plotted the superposition 
of the collective intensity of the BL Lacs and FSRQs for 
the Gilmore and Stecker EBL models and for two dif- 
ferent relative abundances of BL Lacs and FSRQs. For 
the EGRET EGRB data sets, either of the scenarios for 
the relative abundances plotted could reproduce well the 
shapeP^. For the case in which BL Lacs comprise a lower 
fraction of the blazar population, the absorption feature 
is more prominent, as expected given the previous result 
for the LDDE model (Figure O. 

It should be noted that the effects of secondary cas- 
cades have been neglected. Even though the VHE elec- 
trons that are produced in the pair production processes 
of VHE photons can certainly produce VHE photons. 

In principle, in separating the blazar population into its spec- 
tral sub-populations, the parameters for the GLF models would 
have to be refitted to the measured distributions of their corre- 
sponding sub-populations. However, for the sake of discussion, we 
retain the parameters for each GLF model as determined by NT06. 

However, note that there are additional systematic uncertain- 
ties (not plotted here) that could be quite large, especially at the 
highest energies. 
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Fig. 4. — The blazar contribution to the EGRB including EM 
cascade radiation assuming separate emission properties for BL 
Lacs and FSRQs. Thick lines: BL Lac contribution assuming 
the LDDE model and their Fermi-determined SID, (1.97, 0.175). 
Thin lines: FSRQ contribution assuming the PLE model and their 
Fermi-determined SID, (2.37, 0.14). Data points, grey band, solid 
lines, dashed lines, and dot-dashed lines: Same as in Figure [2] 
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Fig. 5. — The superposition of the contributions to the EGRB 
of FSRQs and BL Lacs including EM cascade radiation assuming 
separate GLFs and SIDs for the subpopulations. Solid: The con- 
tribution assuming FSRQs comprise 80% of the blazar population 
and BL Lacs comprise 20%. Dashed: The contribution assuming 
equal numbers of FSRQs and BL Lacs. Data points and grey band: 
Same as in Figure [2] 



the spectra of cascade photons above TeV energies are 
very steep. Since most of the cascade radiation origi- 
nated from the interactions of TeV photons with the in- 
frared part of the EBL, the steep dechne of such primary 
photons in the cascade spectrum ensures that the effect 
of secondary cascades is smaU. We have also neglected 
the possible contribution from synchrotron radiation re- 
sulting from the propagation of VHE electrons through 
intervening magnetic fields, which we will examine in a 
future publication. 

5. DISCUSSION AND CONCLUSIONS 

We have demonstrated the effects of EM cascade radi- 
ation on the contribution to the EGRB attributable to 
blazar spectra finding that the cascade radiation flattens 
the spectrum of the overall collective intensity of blazars 
and can greatly enhance their resulting contribution to 
the EGRB at the high-end of the Fermi energy range. 
We have also demonstrated that the amount of cascade 
radiation is sensitive to the EBL model and the nature of 
blazar spectra. We make use of a distribution of blazar 
gamma-ray spectral indices as determined from the most 
recent data provided by Fermi taking into account the 
possible spread in the distribution. Such a spread in the 



SID results in curvature in the collective intensity, which 
impacts the number of VHE photons available to create 
cascades. Furthermore, we have applied a similar analy- 
sis to two blazar sub-populations, BL Lacs and FSRQs, 
distinguishing between their spectral characteristics and 
their evolution with cosmic time. We found that while 
the inclusion of cascade radiation does not alter much 
the spectrum of the collective intensity of FSRQs, it can 
considerably enhance the collective intensity of BL Lacs. 
Hence, the FSRQs will dominate the collective intensity 
at lower energies and the BL Lacs at higher energies. 
In comparison to the aforementioned analyses (see Sec- 
tion [ij that include semiempirical determinations of the 
cascade radiation, the spectrum of this analysis appears 
different owing to the inclusion of the effect of the spread 
in the blazar spectral index and the distinct considera- 
tions of the evolution and spectral characteristics of the 
BL Lacs and FSRQs. In fact, the distinguishing features 
of this analysis allow the unique outcome of producing a 
collective blazar spectrum that reproduces well the shape 
of the spectrum of the EGRET EGRB. 

As demonstrated in this paper, EM cascades can 
greatly enhance the contribution to the EGRB at- 
tributable to blazars, particularly if the spread in the 
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blazar SID is large or consists of relatively hard mean 
spectral indices and blazar spectra do not change appre- 
ciably at VHEs. As the cascade radiation, as well as the 
emission from the intrinsic spectra, is subject to absorp- 
tion, the overall collective blazar intensity will exhibit an 
absorption feature at the highest energies, the shape of 
which is sensitive to the blazar GLF and the model of 
the EBL. The enhancement of the high-energy intensity 
by cascade radiation allows the high-energy absorption 
feature to be more easily measured by Fermi. On the 
other hand, in the case that BL Lacs and FSRQs form 
distinct blazar sub-populations with respect to spectral 
characteristics and evolution (as indicated bv lAbdo et alj 
l2009al) . the cascade radiation from the BL Lacs could ac- 
count for a significant fraction of the BL Lac contribution 
to the collective blazar intensity, particularly at high en- 
ergies. If the collective intensity from BL Lacs is concen- 
trated at low redshifts, then even though their cascade 
radiation may be considerable, their collective spectrum 
above 10 GeV would exhibit little attenuation. Hence, a 
large fraction of BL Lacs in the blazar population as com- 
pared to the FSRQs could result in a collective spectrum 
for which the absorption feature is less prominent. Thus, 
the observation of the absorption feature by Fermi would 
provide a unique probe into the evolution of blazars with 
cosmic time as well as the nature of the EBL. In fact, as 
shown in the figures, the preliminary Fermi EGRH^ is 
intriguingly steep. However, as this is only a preliminary 
measurement (and without luminosity functions based 
on Fermi blazars and having calculated for the EGRET 
sensitivity), it is, as yet, premature to draw conclusions 
from comparisons with this analysis. Nonetheless, we 
will return to this question in the near future. 

As noted in Section 13.21 it is possible that over the 
course of so many decades in energy, blazar spectra could 
change, most likely turning over at some energy. In such 
a case, the cascade radiation could be reduced consid- 
erably and the collective intensity of blazars might not 
be appreciably enhanced at high energies. The most no- 
table effect is that on the contribution from BL Lacs. 
The flattening of the spectrum at high energies is due, 
at least in part, to the enhanced cascade radiation. If 
BL Lac spectra break at energies low enough to render 
the cascade radiation insubstantial, the collective spec- 
trum at high energies would be steeper. If the break 
occurs at energies greater than 100 GeV, then the col- 
lective spectrum at high energies would still be more en- 
hanced than it would be for the collective intensity from 
FSRQs alone and the absorption feature could be ob- 
servable. If the break occurs at energies on the order 
of tens of GeV, then the collective spectrum would be 
steeper than in the either of the previous scenarios due 
to both the steeping of the intrinsic spectrum and ab- 
sorption. If even without the enhancement from cascade 
radiation blazars comprise the bulk of the EGRB, the 
low-energy-break and high-energy-break scenarios would 
be distinguishable from each other, as well as distinguish- 
able from the no break scenario. Thus, the study of the 
collective intensity of blazars and the measurement of the 
EGRB by Fermi will have implications for the nature of 
blazar spectra at VHEs, which would be complimentary 
to the constraints provided by measurements of individ- 
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ual blazar spectra by lACTs (jAlbert et al.ll2008[ ). 

Additionally, it is possible that blazar spectra change 
not only at high energies, but also du ring periods of 
increased activity ( see, e . g., von Monti gnv et al.l 119951 : 
'Steck er fc SalamonI Il996t iBottcheii Il999f ): hence, the 
SIDs of such flaring blazars would be distinct from that of 
quiescent blazars, affecting their respective contributions 
to the collective blazar intensity. However, the question 
of whether blazar spectra change during flaring periods 
remains, as yet, unresolved due in large part to the lack 
of tim e resolution in EGRET data (Venters & Pavlido3 
|2007| ). With the improved sensitivity of Fermi, time- 
resolved spectroscopy will be possible allowing the mea- 
surement of blazar spectra during flaring periods and the 
determination of whether flaring and quiescent blazars do 
indeed comprise spectrally distinct populations. 

It should be noted that the question of blazar statistics 
and their role in determining the blazar contribution to 
the EGRB is complicated by the presence of interven- 
ing magnetic fields. With negligible magnetic fields, the 
relativistic beaming of the intrinsic gamma-ray emission 
from AGNs ensures that observable gamma-rays origi- 
nate only from those AGNs that are favorably aligned 
with respect to the observer (i.e., blazars). This is true 
both for the intrinsic gamma-ray emission from AGNs 
and also for the cascade radiation, since without mag- 
netic fields, the cascades are highly collimated in the di- 
rection of the propagating VHE photons. In the pres- 
ence of substantial magnetic fields, the cascades will no 
longer be collimated, increasing the number of gamma- 
ray emitting AGNs contributing to the EGRB. The lack 
of collimation of the cascades would also introduce ha- 
los of gamma-ray emission around blazars, which would 
contribute to source confusion. On the other hand, the 
presence of substantial magnetic fields would introduce 
a synchrotron component in the cascades redistributing 
the energy in the spectrum of cascade radiation. Though, 
the study of the effects of intervening magnetic fields 
would be insightful, their inclusion would further compli- 
cate this analysis. Thus, the consideration of magnetic 
fields is out of the scope of this paper; however, we will 
return to this question in a future publication. 

It is also possible that blazars cannot entirely account 
for the EGRB or could account for the EGRB only at 
lower energies. In such a case, the participation of other 
sources of high-energy gamma rays would be required to 
fully explain the EGRB, and the spectrum of the EGRB 
would also exhibit interesting features. For example, 
though blazars currently dominate the identifications of 
extragalactic TeV gamma-ray emitters, it is by no means 
clear that blazars are the only sources of TeV gamma 
rays. If other astrophysical objects, such as GRBs and 
starburst galaxies, are found to also emit TeV gamma 
rays, they too will produce radiation from cascades. 
However, as likely candidate extragalactic TeV sources 
also constitute cosmological populations, they too will 
exhibit absorption at the highest energies. In addition 
to cascades from VHE photons, secondary particles pro- 
duced through the propagation of ultra-high-energy cos- 
mic rays through the EBL and the CMB will initi ate cas- 
cades which will also contribute to the EGRB (iSteckeii 
[19701 IDermedll986tlLei[T998tlStecker fc Salamonll 19991) 7 
The origins of ultra-high-energy cosmic rays are, as yet, 
unknown; hence, the study of the EGRB together with 
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Fermi measurements can have profound imphcations for 
the highest energy cosmic rays. 

On the other hand, if the collective intensity from some 
other gamma-ray source peaks at high energies, then the 
resulting EGRB would also flatten due to the transition 
and could turn over depending on the spectrum of the 
overtaking gamma-ray source. In such a case, Fermi- 
measured blazar GLFs and SIDs and information about 
the high-energy properties of other gamma-ray sources 
will be necessary to distinguish between a transition and 
a higher energy break in blazar spectra. An intriguing ex- 
ample to consider would be that of gamma-ray emission 
from dark matter annihilation, which would peak near 
the mass of the dark matter particle. If dark matter an- 
nihilation contributes substantially to the EGRB (even if 
only at energies near the mass of the dark matter parti- 
cle), then i t will be revealed t hrou gh anisotropy searches 
(see, e.g., lAndo et all l2007t iSieeal-Gaskins &: Pavlidoui 



MM)- 

Thus, the study of the collective blazar intensity, in- 
cluding the effects of cascade radiation, combined with 
the measurements of the EGRB and the blazar gamma- 
ray luminosity-redshift distribution by Fermi can provide 
a wealth of insight into the nature of blazar spectra, as 
well as the evolution of blazars and the EBL with cosmic 
time. 
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APPENDIX 

DERIVATION OF THE SPECTRUM OF ELECTROMAGNETIC CASCADE RADIATION 
To determine the spectrum of cascade photons, we begin with the spectrum of scattered radiation per electron: 



dNr 
'dE' 



dtdEde 



-dtde 



r' 



dtdEde 



dt 



dE, 



dEede , 



(Al) 



where F is the initial Lorentz factor of the electron, F' is the Lorentz factor of the electron at some later time t, 
and we are integrating over the soft photon background and the inverse Compton timescale, which is determined by 
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the electron energy losse^Fl. \dt/dEe\. The differential scattered p hoton spectrum is determined by integrating the 
differential Klein-Nishina cross section (jBlumenthal k, Gouldlll970[ ): 



d^NT' Sarc 7i(e) 



dtdEde 4r'2 e 



1 (r.g) 



2 



2qlnq+{l + 2q) {1 - q) + 3 " 



(A2) 



where n(e)de is the soft photon background number densit\F'1. = ieT'/mc^, 

q = rj/T^{l — rj), 77 = ei/T'mc^, and ei is the scattered photon energy. The spectrum of scattered photons 
also defines the electron loss rate: 

since adiabatic losses due to cosmic expansion that occur between interaction points for the electron are negligible 
compared to the inverse Compton losses. The spectrum of inverse Compton scattered photons arising from the 
distribution of electrons per primary photon is determined by integrating over the spectrum of electrons, dN^ = 
NeiT)dr: 



dE J dE J dE 

We can determine the cascade spectrum by integrating over the spectrum of primary photons, dN^ — Ny{Ep)dEp 

dNtot 



dE 



(A5) 



However, in the case of any photon interaction process that cannot be considered continuous, the "spectrum" of 
electrons per primary photon is discrete. In this case, the resulting spectrum of inverse Compton scattered photons is 
best described as 

dE-^ dE ■ ^^^> 

For a discrete distribution of primary photons (i.e., Monte Carlo photon propagation), the total cascade spectrum is 
given by 

N'' 



where A^ is the number of electron-like particles produced in the interaction of the k*"^ primary photon, and Tj is 

the Lorentz factor of the j*'^ electron resulting from the interaction of the fc**^ primary photon. In the case of pair 
production, N^^ = 2, Fj = / x Ep./meC^, and F^ = (1 — /) x Ep^/nieC^, where / is some fraction of the energy of 
the primary photon (determined from the kinematics of the interaction). The observed spectrum from a particular 
epoclH would then be 



dEo ^ ' ^ ' ^ dE 



= + 2. l.^- (A8) 
fc=i i=i 



The total observed spectrum is determined by summing over all of the epochs: 

" 1=1 fe=i j=i 

where A'zbins is the number of redshift bins, zi is the redshift of the l^^ bin, and Az; is the size of the Z*^ redshift 
bin. For a continuous spectrum of photons undergoing pair production interactions, d'^N^/dzdEp, the total observed 
spectrum is given by 



-(Eo)^ + z)—-—F{},Ep,z) 



dEo J J dzdEp 



dE dE 



dEpdz , (AlO) 



At the energies that we are considering in this paper, the continuous energy loss description for electrons is accurate even at energies 
that are no longer quite in the Thomson regime. However, in the extreme Klein-Nishina regime, electrons can lose large fractions of their 
energies in individual interactions rendering the continuous energy loss approximation invalid. 

Note that the soft photon background is redshift dependent. However, the inverse Compton timescale is sufficiently short so that the 
soft photon background does not change appreciably during the creation of the cascade. 

Unlike electrons, photons can propagate a considerable distance before the point of interaction and adiabatic losses must be included. 
Hence, Ep is the redshifted energy of the primary photon. 
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where P(/; Ep, z) is the probability that the pair production interaction of a primary photon of energy Ep will produce 
electron-type particles of energies / x Ep and (1 — /) x Ep. Including the attenuation factor, exp[— t(£^Oj z)], we finally 
arrive at the observed diiferential cascade spectrum: 



'dEo 



dNr,{Eo{l + z)) 



d'^N. 



dE 

dN,jEMl + z)) 



dE 



e-^{Eo,^)dE^dz 



(All) 



